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Abstract 

In this paper, we consider the two-loop sunset diagram with two different masses, 
m and M, at spacelike virtuality q 2 = —m 2 . We find explicit representations for 
the master integrals and an analytic result through 0{e) in d = 4 — 2e space-time 
dimensions for the case of equal masses, m = M. 
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1 Introduction 



The two-loop self-energy diagram with three lines, the so-called sunset diagram, plays a 
very important role in the evaluation of higher-order corrections in quantum electrody- 
namics (QED), quantum chromodynamics (QCD), and the electroweak theory. In fact, it 
enters almost every calculation of such kind. Therefore, it has been an object of intensive 
investigations in the literature (see, e.g., Ref. jT]). In the cases when one or two masses 
are zero, analytical results in terms of polylogarithmic functions are available [21 0] even 
at arbitrary external momentum q. 

A much more difficult situation arises when all three lines are massive. In this case, the 
result cannot be expressed in terms of polylogarithms, and more involved functions must 
be considered j3]. In particular, the case when all three masses are equal, m\ = m-i = 
m3 = m, is of great importance. Calculations for such kinematics were carried out using 
the dispersion-relation representation [Sj, the differential-equation approach, and several 
expansions in the variable q 2 /m 2 [H]- Recently, a detailed analysis of the differential- 
equation approach, together with solutions and analytical continuations, was presented 
[7j. At singular points, including zero invariant mass (q 2 = 0), threshold (q 2 = 9m 2 ), and 
pseudothreshold (q 2 = m 2 ), the values of the diagram can be found in closed form. At 
any other value of q 2 , we have formal integral representations. There is, however, another 
interesting point where one would like to have some closed analytical expression, namely 
at q 2 = —m 2 . Such kinematics occurs, for example, in the calculations of hard Wilson 
coefficients in non-relativistic QED and QCD (see, e.g., Ref. |B]) at the next-to- leading 
order (NLO). This can be applied, e.g., to evaluation of the parapositronium decay rate 
in non-relativistic QED. In non-relativistic QCD, applications are related to the analysis 
of heavy- quarkonium decays and the near-threshold ti production at a future e + e~ linear 
collider, which offers the opportunity to improve the accuracy of the t-quark mass. 
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Figure 1: (a) Two-loop sunset diagram J^^ M and (b) one- loop self-energy diagram I^ M 
with different masses, m and M. In both cases, the external four-momentum q satisfies 
the condition q 2 = —m 2 . The label attached to a given line indicates the mass square 
appearing in the respective propagator and the power to which the latter is raised. 

The goal of this paper is to fill this gap by finding a suitable representation for the 
above-mentioned object. Actually, we solve an even more general problem. In fact, we 
evaluate the sunset diagram with two different masses, m and M, at q 2 = —m 2 , as shown 
in Fig.^a). The case of interest then simply emerges by setting m = M. 
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This paper is organized as follows. In Section we establish a relation between the 
two- loop sunset diagram of Fig. [TJa) and the one- loop one of Fig. H£b). In Sectional we 
work out the one-loop case. In Section 01 we consider the two-loop case for the general 
case of m ^ M. In Section El we derive a closed expression for the special case of m = M. 
Our conclusions are contained in Sectional 



2 Relation between one- and two-loop sunset dia- 
grams 

In our calculation, we use a particular formula to represent a loop with two massive 
propagators as an integral whose integrand contains a new propagator with a mass that 
depends on the variable of integration. Graphically, this formula has the following form: 



a,M 2 




■l+d 



b,AI 



T(a + b-d/2) 

r(a) r(6) 



ds 



'\ — s ")a+l-<2/2 s &+l-d/2 



a+b-d/2 . . 

X „ * „ * , 1 

Q M 2 /[s(l-s)] ' V ' 

where the loop with two propagators with mass square M 2 is replaced by one propagator 
with mass square M 2 /[s(l — s)]. Here and in the following, we employ dimensional 
regularization in d = 4 — 2e space-time dimensions. 

Equation (JTJ) is easily derived from the Feynman parameter representation and was 
introduced in Euclidean and Minkowski spaces in Refs. and [TU], respectively. Here, 
we work in Minkowski space and will, thus, follow Ref. |l(Jj . 

We shall adopt the following strategy: applying Eq. (JTJ), we shall represent the results 
for the sunset diagrams as integrals of one-loop two-point diagrams involving propagators 
with masses that depend on the variable of integration. We shall first evaluate the one-loop 
integrals in Section El and then reconstruct the results for the two-loop sunset diagrams 
involving two different masses with the help of Eq. in Section A similar strategy 
was recently also adopted for the calculation of certain four-loop tadpole diagrams |TT] . 



3 One-loop case 

Let us first consider the one-loop case shown in Fig. ^b). The self-energy integral in d 
dimensions is given by 

jmM _ f d^fc 1 

ab J TT d / 2 (k 2 - m 2 ) a [{k - q) 2 - M 2 } b ' 1 ' 

where q 2 = —m 2 is implied. 
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First we try to establish the general structure of the result. For this purpose it is 
convenient to consider the integral I^ M rather than I™ . This integral is ultraviolet 
finite, so that we may put d = A. (We will return to the case d = 4 — 2e later on.) Using 
the standard Feynman parameter technique, it is not difficult to obtain the representation 



jmM 
'12 



—I 



dt (1 - t) 



t)M 2 + t(2 - t)m 2 



(3) 



which could be further evaluated in terms of logarithms. 

Having in mind the algorithm proposed above, our next step will be to replace the 
mass M 2 by M 2 /[s(l — s)] with s integrated from to 1. The easiest way to perform such 
an integration is to expand Eq. Q in powers of M 2 or 1/M 2 . Then, the integration over 
s can be performed term by term leading to Euler gamma functions. We shall follow this 
logic and expand I^ M in 1/M 2 . The expansion in M 2 , corresponding to another regime, 
is also possible, but more involved and will not be considered here. 

It turns out that Eq. (jSJ) is not appropriate for obtaining the coefficients of the ex- 
pansion. It is much more advantageous to start from a differential equation. First, a 
differential equation for I™ M is obtained through integration by parts (for a review of this 
procedure, see Refs. j2JH2])- Next, differentiating both sides of this equation w.r.t. M 2 , 
one obtains a differential equation for I^ M , exploiting that I^ M = (d/dM 2 )/{^ M . This 
equation reads 



d-A + 



Am 4 



1 + 



4m 4 
M 4 " 



dM 2 



jmM 
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12 



(4) 



where jyi is a non-uniform term which contains only vacuum tadpoles and is given by 



f ( m \ 1 

/l2 [M 2 J ~ 



m 



}f 2 'M i \ lU M 2+1 , 



+ 0(e) 



(5) 



In order to solve the differential equation (|4*jl. we need one boundary condition. A suitable 
condition follows from the observation that I^ M — > as M 2 — > oo. Then, the solution of 
Eq. flU) reads 



m j 



mM i 
12 



dx' m 2 

—fi2{x) 



VI + 4x 2 J x'y/l + Ax' 2 i 



dx' 



VTTIx 2 J VTTax 12 



'l + 2x' + 2x'lnx') 



(6) 



where d — 4 and we have introduced the variable 

-,2 



X 



m 
M 2 



(7) 
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The first term in Eq. (JUJ) can be rewritten as 

dx' —x „ /l 1 3 



F,[ -Ax 2 



VI + 4x 2 J y/1 + Ax' 2 y/1 + 4x 2 V2'2'2 

where 2 -Pi is the hypergeometric function Using the property of the hypergeometric 
function [T5] 

2 F 1 ((i,6;c;x) = (l-i) c - fl - fc 2 F 1 (c-a,c-5;c;i), (9) 
Eq. JSJ can be transformed into 

-^(1,1;-; -4* 2 ) = - £ (-16x 2 )^ - . (10) 



n=l 



The sum of the last two terms on the r.h.s. of Eq. © yields 



{l ' Jj " - ; (\nx' + l) = U - 1 -l)lnx (11) 



VI + 4x 2 J VI + 4a;' 2 2 V VI + 4a; 2 



2\n 



g r(n + l/2) (-4a; 2 ) 



4 VI + 4x 2 n!r(l/2) n 
The term in the second line of Eq. (jllj) can be represented as a limit, 



lim - 

<5^0 5 



2 Aj (^,5; 1 + 5; -4a; 2 ) - I 



(12) 



Now, using property (jHJ), we can remove the square root VT + 4a; 2 . Finally, after some 
transformations, we obtain the following series representation for the one-loop diagram 

mil 
12 



J™ M with d = 4: 



=V<») - £ :(-*")■ ( 2 D (5 m» - \s l{n - 1) + i Sl(2n - i) - JL 



+ L? ( - ifa2) "^)^' (13) 



n=l 



where S'i(n) = Sj=i 1/j is the harmonic sum. This is just the form we wanted to have. 
The appearance of the central binomial coefficient here is not surprising. It corresponds 
to the branch point at x — 4. Similar one- and two-loop diagrams were considered in 
Ref. [H], and binomial sums were studied in Ref. |To] . 

Now, with the knowledge of the generic n-th coefficient of the expansion in x, it is not 
difficult to find subsequent terms in the e expansion. To do this, one has to increase the 
"weights" of the harmonic sums arising from 1/n, S± to 1/n 2 , Sf, and S2 and so on. At 
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0(e), we thus obtain 



m T mM,0(e) 
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(14) 



4 Two-loop sunset diagram 

Now, we are ready to evaluate the two- loop sunset diagram J^ IM shown in Fig. Q^a). 
Its loop-integral representation reads 

jmMM _ f d d l 1 /-fK\ 

abc "7 7r d (P-m 2 ) a (Z 2 -M 2 ) 6 [(/-£;-g) 2 -M 2 ] c ' 1 J 

where g 2 = —m 2 is implied. 

Exploiting again Eq. we replace M 2 in the one-loop subdiagram by M 2 /[s(l — s)] 
and integrate over s. It is known, e.g. from Ref. ^H]; that for such kinematics there exist 
three master-integrals. As master integrals we choose J[^^ M , J^^ M , and J]T^ M . 

After some calculation, we finally obtain 

M 2jmMM = 1 yv 2\n ( n ) / ]_ j #1 _ 3 5j. Si 1 

z±T W 4 «M 2n 2n 2n n An 2 



71=1 



2;; 



1 ~ (-16a: 2 )" / 1 J_ 1 \ 
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x£±x ' (f) n(2n- l) 2 V2 2 2 

1 1 1 \ . « (-16x 2 )" -1 

^ f 2 ») Q") 4n 2 (2n+l) ' 
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+ x V(-x 2 )"ifi i , ( - Inx - -S 1 + -ST - S~i 
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1 1 1 \ 



+ E 



4n 4(n + 1) An + 1 y 
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n=l 



(IKS) 2n 2 (2n + l)(2n 



1) ' 



where ( denotes Riemann's zeta function and, for brevity, we omitted the arguments of 
the harmonic sum Si(n) and introduced the short-hand notations 

S 1 = S 1 {n-l), = Si(2n - 1) , ^=5i(4n-l). (19) 

For the practical applications mentioned in the Introduction, we also need the 0(e) 
terms of the sunset diagrams. We do not present them here for a general value of x, since 
the corresponding expressions are cumbersome. However, we shall do this in the next 
sections for the equal- mass case, which is of most interest. 



5 Equal- mass case 

We now turn to the case when m = M. Note that, with equal masses, there are only two 
master integrals, J m and Jn 2 . After substituting x = 1 in Eqs. (|16 p -()18j ) . one obtains 
series which cannot be summed explicitly to become some known constants. However, 
they can be slightly simplified. In fact, using the PSLQ algorithm J7|, which is able to 
reconstruct the rational- number coefficients multiplying a given set of irrational numbers 
from a high-precision numerical result, we find that there are some relations among sums 
at x = 1 and the irrationals In 2 and C(2)- Specifically, we can establish the following. 
Among the sums occurring in Eqs. (fTB |) -([T% |) . 



/2n\ 
]\n V n J 



1111 1 1 1 



| ^ [n* n 2 ' ri + 1' (n + 1) 2 ' 2n-l' (2n-l) 2 ' 4n-l' 

1 1 (p (p 4> cf> 1 

An + 1' (An + l) 2 ' n' n + l' 2n- l' (2n - l) 2 ' 4n + 1 J ' 

where = 5"i — 3Si + 2Si, are not independent. To this list, we add also the corresponding 
sums with the prefactors 1 or 0/n 2 . All of them can be expressed in terms of five sums 
and the constants In 2 and C(2)- We are able to evaluate two of these sums in terms of 
known functions. Let us introduce the notation 

p — Vb } (p = arctan y/p . (20) 
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Then, the new constants read 



fx = ^pF [2 arctan y/p, ^— ^ I = —pF [ 2cp, — V~ ) = 1.8829167613 . . . , (21) 



1+p^ 




2p y 








2p ; 





d = — £ ( 2 arctan ) = — E ( 2<p, -=— ) = 0.9671227369 . . . , (22) 

v/P V 2p J VP \ 2sin>/ 

where F and E 1 are the elliptic functions of the first and second kind, respectively. For 
the reader's convenience, we recall here their definitions: 

F(<f>,k) = / d0(l-A;sin 2 0)- 1/2 , E{6, k) = [ d6 (1 - k sin 2 6) 1/2 . (23) 
Jo Jo 

With these constants, we evaluate the following two sums: 

( 2n ) 111 1 

£ (-irkr = -2 + To p + s pei - To p/l ' (24) 



(an) 



^ ^->ff = — 1H2 + -A. (25) 

n (S) 3 3 



n=l 



It seems that other possible irrationalities will include some integrals of elliptic functions. 
This requires further investigations. 

Among the five sums of the other kind appearing in Eqs. (|16| ) - (fT8|h 

y (-16)" f 1 l_ 1 1 ) ^ 

k( 2 :)(Z)r n (2n-lY (2n + l)j > [ b) 

only two are independent. Here is an example of a relation between these sums: 



n=l 



( 2 :) ( 



2n 



It is interesting to observe that sums of different "weights," namely 1, 1/n, and 1/n 2 , 
appear here. In earlier studies, the independence of sums with different "weights" was 
encountered and frequently used (see, e.g., Refs. [9*1 ITU1 ITH]). 

A similar analysis can also be performed at 0(e). Finally, setting m = M = 1, we 
obtain for the master integrals through 0(e) 

Sm + eE^ + Ofc 2 ), (28) 
J™ = ~ - 77; + S 112 + eY$$ + 0(e 2 ) , (29) 
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where 

215 9 , x 

= 24" + 4 ' ' 

'2n 



S? 1 ( r ) = C(3) + yC(2)-yC(2)ln2 ]jj 



-9.03056576107922587907436223954936770213033473413 . . . , (30) 



Sll2 = ^ + i c(2) 

18 4 1 ; 



0.0113804720812563731826135489564394881100859890024139 
49 w _ 15 „ 1607 



2» 



+ 



^ 65 5 5 , 77 15 o 15 0p 77 15 k 

S (£) V T "8 (7 "4 #+ Tn + T6n + Y- + T6^ + T^ 

45 \ ~ (-16) n /655 50 493 25 /3\ 



1G/i 3 J + 5 M f 4n "l V 18 + 3 13 3C// 3 /; 



n=l 



\ n J \2n 



= -25.4473853869869254012665862897193562434443269569232788372 . . . , (31) 
^ = ic(3) + >)-^.n 2 + f 



n^l ^ ft) l 2 



8 4 r n 16 n 8 n 4 n 2n 2 



+ 



2n, 

15 \ ~ (-16) n /470 100 107 35 /5 s 



/2n\ /< 

= -2.04235726589417810442922843329658331200829586895 .... (32) 

Here, we have introduced the following short-hand notations 

= S 1 {n - 1) - 3Si(2n - 1) + 2S'i(4ra - 1) , 
p = - 1) + 55i(2n - 1) - 2Si(4n - 1) , 
a = S 2 (n - 1) + 10S 2 (2n - 1) - 8S 2 (4n - 1) , 
K = S^n - 1) - 5 1 (2n - 1) + Si(4n - 1) , 

/3 = 2Si(2n-l)-Si(4n-l). (33) 

6 Conclusions 

In this paper, we considered the two-loop sunset diagram with two different masses, m and 
M, in the special kinematical regime where q 2 = —m 2 . The coefficients of the expansion 
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in the variable x = m 2 /M 2 were found explicitly. This allows for the numerical restoration 
of the result, for example via Pade approximants. In the special case of m = M, the result 
was expressed in terms of very fast-converging alternating series. Some analytical relations 
between these series were found using the PSLQ algorithm P^J- Their exact analytical 
structure is related to elliptic integrals. The two lowest elements of this structure were 
expressed in terms of the elliptic functions E and F. Integral representations for the sunset 
master integrals J™f IM , J^l^ M , and J^^ M will be presented in a future publication. 

Acknowledgments 

The authors are grateful to M.Yu. Kalmykov for fruitful discussions and valuable 
comments. The work of A.V.K. and O.V. was supported in part by RFBR Grant No. 
N 05-02-17645. The work of A.V.K. was also supported in part by the Heisenberg-Landau- 
Programm. This work was supported in part by BMBF Grant No. 05 HT4GUA/4 and 
HGF Grant No. NG-VH-008. 



References 

[1] P. Post, J.B. Tausk, Mod. Phys. Lett. A 11 (1996) 2115; 

FA. Berends, A.I. Davydychev, N.I. Ussyukina, Phys. Lett. B 426 (1998) 95; 

S. Groote, J.G. Korner, AA. Pivovarov, Phys. Lett. B 443 (1998) 269; 

S. Groote, J.G. Korner, AA. Pivovarov, Nucl. Phys. B 542 (1999) 515; 

S. Groote, J.G. Korner, AA. Pivovarov, Eur. Phys. J. C 11 (1999) 279; 

A.I. Davydychev, VA. Smirnov, Nucl. Phys. B 554 (1999) 391; 

A. Bashir, R. Delbourgo, M.L. Roberts, J. Math. Phys. 42 (2001) 5553; 

P. Mastrolia, E. Remiddi, Nucl. Phys. B 657 (2003) 397; 

A. Onishchenko, O. Veretin, Yad. Fiz. 68 (2005) 1461 [Phys. Atom. Nucl. 68 (2005) 
1405]. 

[2] A.V. Kotikov, Phys. Lett. B 254 (1991) 158; 
A.V. Kotikov, Mod. Phys. Lett. A 6 (1991) 677. 

[3] J. Fleischer, F. Jegerlehner, O.V. Tarasov, O.L. Veretin, Nucl. Phys. B 539 (1999) 
671; 

J. Fleischer, F. Jegerlehner, O.V. Tarasov, O.L. Veretin, Nucl. Phys. B 571 (2000) 
511, Erratum; 

A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B 605 (2001) 266; 

F. Jegerlehner, M.Yu. Kalmykov, O. Veretin, Nucl. Phys. B 658 (2003) 49. 

[4] FA. Berends, M. Buza, M. Bdhm, R. Scharf, Z. Phys. C 63 (1994) 227; 

S. Bauberger, FA. Berends, M. Bdhm, M. Buza, Nucl. Phys. B 434 (1995) 383. 

[5] D.J. Broadhurst, Z. Phys. C 47 (1990) 115; 

A.I. Davydychev, R. Delbourgo, J. Phys. A 37 (2004) 4871. 



10 



[6] D.J. Broadhurst, J. Fleischer, O.V. Tarasov, Z. Phys. C 60 (1993) 287. 



[7] S. Laporta, E. Remiddi, Nucl. Phys. B 704 (2005) 349. 

[8] A. Czarnecki, K. Melnikov, A. Yelkhovsky, Phys. Rev. A 61 (2000) 052502; 
A. Czarnecki, K. Melnikov, Phys. Rev. D 65 (2002) 051501. 

[9] J. Fleischer, A.V. Kotikov, O.L. Veretin, Phys. Lett. B 417 (1998) 163. 

[10] J. Fleischer, A.V. Kotikov, O.L. Veretin, Nucl. Phys. B 547 (1999) 343. 

[11] B.A. Kniehl, A.V. Kotikov, Report No. DESY 05-153, |hep^ph/0508238| 

[12] A.V. Kotikov, Phys. Lett. B 259 (1991) 314; 
A.V. Kotikov, Phys. Lett. B 267 (1991) 123; 

E. Remiddi, Nuovo Cim. A 110 (1997) 1435. 

[13] M. Abramowitz, LA. Stegun (Eds.), Handbook of Mathematical Functions with For- 
mulas, Graphs, and Mathematical Tables, National Bureau of Standards Applied 
Mathematics Series, vol. 55, US Government Printing Office, Washington, DC, 1964. 

[14] F. Jegerlehner, MYu. Kalmykov, O. Veretin, Nucl. Phys. B 641 (2002) 285; 

F. Jegerlehner, MYu. Kalmykov, O. Veretin, Nucl. Phys. B 658 (2003) 49. 

[15] MYu. Kalmykov, O. Veretin, Phys. Lett. B 483 (2000) 315; 

A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B 699 (2004) 3. 

[16] O.V. Tarasov, Nucl. Phys. B 502 (1997) 455. 

[17] H.R.P. Ferguson, D.H. Bailey, RNR Technical Report No. RNR-91-032; 

H.R.P. Ferguson, D.H. Bailey, S. Arno, NASA Technical Report No. NAS-96-005. 

[18] J. Fleischer, M.Yu. Kalmykov, A.V. Kotikov, Phys. Lett. B 462 (1999) 169; 
A.V. Kotikov, L.N. Lipatov, Nucl. Phys. B 582 (2000) 19; 
A.V. Kotikov, L.N. Lipatov, Nucl. Phys. B 661 (2003) 19; 

A.V. Kotikov, L.N. Lipatov, A.I. Onishchenko, V.N. Velizhanin, Phys. Lett. B 595 
(2004) 521; 

Z. Bern, L.J. Dixon, D.A. Kosower, Comptes Rendus Physique 5 (2004) 955; 
Z. Bern, L.J. Dixon, V.A. Smirnov, Phys. Rev. D 72 (2005) 085001; 
S. Moch, J.A.M. Vermaseren, A. Vogt, JHEP 0508 (2005) 049; 
T. Gehrmann, T. Huber, D. Maitre, Phys. Lett. B 622 (2005) 295. 



11 



